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Q . scope of a Bianchi type I (BI) gravitational field in presence of a viscous fluid and A term 

CN \ has been carried out. The system of equations defining the evolution of the volume scale 

&JQ' of BI universe, energy density and corresponding Hubble constant has been derived. The 

system in question has been thoroughly studied qualitatively. Corresponding solutions are 



graphically illustrated. The system in question is also studied from the view point of blow 
up. It has been shown that the blow up takes place only in presence of viscosity. 
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The problem of an initial singularity still remains at the center of modern day cosmology. 
Though the Big Bang theory is deeply rooted among the scientists dealing with the cosmology of 
the early Universe, it is natural to reconsider models of a universe free from initial singularities. 
Another problem that the modern day cosmology deals with is the accelerated mode of expansion. 
In order to answer to these questions a number of theories has been proposed by cosmologists. It 
has been shown that the introduction of a nonlinear spinor field or an interacting spinor and scalar 
fields depending on some special choice of nonlinearity can give rise to singularity free solutions 
in one hand [1, 2, 3, 4], on the other hand they may exploited to explain the late time acceleration 
[5, 6]. 

Why study a nonlinear spinor field? It is well known that the nonlinear generalization of clas- 
sical field theory remains one possible way to overcome the difficulties of a theory that considers 
elementary particles as mathematical points. In this approach elementary particles are modelled 
by regular (soli tonlike) solutions of the corresponding nonlinear equations. The gravitational field 
equation is nonlinear by nature and the field itself is universal and unscreenable. These properties 
lead to a definite physical interest in the gravitational field that goes with these matter fields. We 
prefer a spinor field to scalar or electromagnetic fields, as the spinor field is the most sensitive to 
the gravitational field. 

Why study an anisotropic universe? Though spatially homogeneous and isotropic, Friedmann- 
Robertson-Walker (FRW) models are widely considered as a good approximation of the present 
and early stages of the Universe. However, the large scale matter distribution in the observable 
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Universe, largely manifested in the form of discrete structures, does not exhibit a high degree of 
homogeneity. Recent space investigations detect anisotropy in the cosmic microwave background. 
The Cosmic Background Explorer's differential radiometer has detected and measured cosmic 
microwave background anisotropics at different angular scales. 

These anisotropics are supposed to contain in their fold the entire history of cosmic evolution 
dating back to the recombination era and are being considered as indicative of the geometry and 
the content of the Universe. More information about cosmic microwave background anisotropy 
is expected to be uncovered by the investigations of the microwave anisotropy probe. There is 
widespread consensus among cosmologists that cosmic microwave background anisotropics at 
small angular scales are the key to the formation of discrete structures. The theoretical argu- 
ments [7] and recent experimental data that support the existence of an anisotropic phase that 
approaches an isotropic phase leads one to consider universe models with an anisotropic back- 
ground. 

Why study a system with viscous fluid? The investigation of relativistic cosmological mod- 
els usually has the energy momentum tensor of matter generated by a perfect fluid. To consider 
more realistic models one must take into account the viscosity mechanisms, which have already 
attracted the attention of many researchers. Misner [ , 8] suggested that strong dissipative due to 
the neutrino viscosity may considerably reduce the anisotropy of the black-body radiation. Viscos- 
ity mechanism in cosmology can explain the anomalously high entropy per baryon in the present 
universe [9, 10]. Bulk viscosity associated with the grand-unified-theory phase transition [ I ] may 
lead to an inflationary scenario [12, 13, 14]. 

A uniform cosmological model filled with fluid which possesses pressure and second (bulk) 
viscosity was developed by Murphy [15]. The nature of cosmological solutions for homoge- 
neous Bianchi type I (BI) model was investigated by Belinskii and Khalatnikov [16] by taking into 
account dissipative process due to viscosity. They showed that viscosity cannot remove the cos- 
mological singularity but results in a qualitatively new behavior of the solutions near singularity. 
They found the remarkable property that during the time of the big bang matter is created by the 
gravitational field. 

Given the importance of both viscous mechanism and nonlinear spinor field we have recently 
studied the system in question from various aspects. In [17] we have studied the evolution of a 
BI universe filled with viscous fluid in presence of a A term. Exact solutions to the corresponding 
system of equations were found for some special choice of viscosity parameters. This study was 
further developed in [18], where the system was studied qualitatively. Introduction of a nonlin- 
ear spinor field into the system considerably changes the situation giving rise to some unexpected 
results such as Big Rip without phantom dark energy. The system in question was studied ana- 
lytically in [19, 20] and generalized in [21] employing both numerical and qualitative methods. 
Since the interacting system of spinor and scalar fields gives rise to a induced nonlinearity of the 
spinor field that can change the picture drastically, we plan to consider this system as well. Some 
exact solutions to the system of equations were obtained in [22]. Here we thoroughly study the 
interacting spinor and scalar fields within the framework of a BI gravitational field in presence of 
a viscous fluid and A term. In doing so we will exploit both numerical and qualitative methods. 

II. BASIC EQUATIONS 

We consider a self-consistent system of interacting nonlinear spinor and scalar fields within 
the scope of a Bianchi type-I (BI) gravitational field filled with a viscous fluid in presence of a 
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cosmological term. Corresponding Lagrangian takes the form: 



. m ^r VA+ _ ( p ja( p-«(l + XF), 



(2.1) 



Here m is the spinor mass, A is the coupling constant and F = F(I,J) with I = S 2 = (Y¥) 2 m & 
J = P 2 = (ixffy 5 y) 2 • According to the Pauli-Fierz theorem among the five invariants only I and J 
are independent as all other can be expressed by them: Iy = =7+7 and Iq = I — J. Therefore, 
the choice F = F(I,J), describes the nonlinearity in the most general of its form [3]. Note that 
setting A = in (2.1) we come to the case with minimal coupling. 

The gravitational field in our case is given by a Bianchi type I (BI) metric 



ds 2 



2 2 2 

dt —adx 



b 2 dy 2 -c 2 dz 2 , 



(2.2) 



with a, b, c being the functions of time t only. Here the speed of light is taken to be unity. 
For the BI space-time (2.2) on account of the A term this system has the form 

b c be i 

- + - + -- = kT\ + A, 

b c be 



c a c a 

- + -H 

c a c a 

d b db 

a b ab 

db be cd 

ab be ca 



kT 2 + A, 



kTJ + A, 



kT^ + A, 



(2.3a) 
(2.3b) 

(2.3c) 

(2.3d) 



where over dot means differentiation with respect to t and Ty is the energy-momentum tensor of 
the material field given by 



TP — ' nP V 

V - 4 



WYn v v W + WYv V M W - V M yry v y - V v i//y M y/ 
+ {l-XF) W P-8 p ^ + T jX \. 



(2.4) 



Here TX m is the energy-momentum tensor of a viscous fluid having the form 



1 jlm 



(e + p')u ll u v - p'8^ + T]g v P [u^.p + up.^ - u^u a up. a - upu a u^ a ] 



where 



P = P ~ {£, - 



(2.5) 



(2.6) 



Here £ is the energy density, p - pressure, r\ and £, are the coefficients of shear and bulk viscosity, 
respectively. In a comoving system of reference such that = (1, 0, 0, 0) we have 
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(2.7a) 
(2.7b) 

(2.7c) 
(2.7d) 
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WE consider the case when both the spinor and the scalar fields depend on t only. We also define 
a new function 

T = abc, (2.8) 

which is indeed the volume scale of the BI space-time. It was shown in [19, 20, 22] that the 
solutions of the spinor and scalar field equations can be expressed in terms of x. Then for the 
components of the energy-momentum tensor we find 



mS + 



C 2 



@S + &P- 



C 2 



+ £ 



<0 > 



a 



2t 2 (1+AF) 
C 2 

2t 2 (1+AF) 
C 2 

2t 2 (1+AF) 



p ' + 271 - = T{+2t 1 -, 



a 



.y+2T]- = f 1 1 +2T]-, 



y + 2T] - = 7y+2T7-, 



In account of (2.9) from (2.3) we find the metric functions [3] 

e -2K fridt 



a(t) 
bit) 
c(t) 



yiT 1/3 exp 
y 2 T 1/3 exp 
F 3 T 1/3 exp 



3 

T 
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T(0 
e -2Kjrjdt 

e -lKfr\dt 



with the constants Y; and Xj obeying 

YiY 2 Y 3 = l, 



X l +X 2 + X 3 = 0. 



(2.9a) 
(2.9b) 
(2.9c) 
(2.9d) 



(2.10a) 
(2.10b) 
(2.10c) 



As one sees from (2.10a), (2.10b) and (2.10c), for % = t" with n > 1 the exponent tends to unity 
at large t, and the anisotropic model becomes isotropic one. 

So one needs to find the function T, explicitly. Corresponding equation can be derived from 
Einstein equations and Bianchi identity [a detailed description of this procedure can be found in 
[19, 20, 22]]. For convenience, we also define the generalized Hubble constant. The system then 
reads [22]: 



t = 3//T, (2.11a) 
H = |(3 { H_.)_(3 H «-„-A) + f(= + 5 ^), (2,1b) 

i = 3H(3iH - CO) +4n (3H 2 - KB - A) -Ar\ [k(™ + ^ + )]. (2.11c) 

Here k is the Einstein's gravitational constant, A is the cosmological constant, A is the self- 
coupling constant, m is the spinor mass and n is the power of nonlinearity of the spinor field (here 
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we consider only power law nonlinearity). In (2.11) r\ and £, are the bulk and shear viscosity, 
respectively and they are both positively definite, i.e., 

tj >0, £>0. (2.12) 

They may be either constant or function of time or energy. We consider the case when 

r}=Ae a , %=B£ p , (2.13) 
with A and B being some positive quantities. For p we set as in perfect fluid, 

p=££, C<E(0,1]. (2.14) 

Vismpl05 Note that in this case £ ^ 0, since for dust pressure, hence temperature is zero, that 
results in vanishing viscosity. Note that a system in absence of spinor field has been studied in 
[17, 18]. In that case the corresponding system is analogical to the one given in (2.1 1) without the 
third terms in (2. 1 lb) and (2.11c). 

HI. QUALITATIVE ANALYSIS 

The study of the behavior of dynamic system given by a system of ordinary differential equa- 
tions implies the survey of all possible scenarios of development for different values of the problem 
parameters. It is necessary to understand at least how the process of evolution comes to an end if it 
does so at infinitively large time for a given set of initial conditions which can be given anywhere. 

So, under the specific behavior of the system we understand the phase portrait of the system, 
i.e., the family of integral curves, covering the total phase space. It is easy to imagine as far as 
any point of the space can be declared as the initial one and at least one integral curve will pass 
through it (or it will be fixed point). 

Certainly, it is difficult to imagine such a set of curves. In many cases, close (and not only) 
curves transform into each other at some diffeomorphism of space. These curves are known as 
topologically equivalent. The differences between them are not very important for our study. 
They all behave in the same manner. This relation - "the relation of equivalence" - divides the 
family of curves into the classes of equivalence. For graphical demonstration it will be convenient 
to present at least one representative of each class. 

The change of the value of problem parameters not always results in significant change of the 
phase portrait. Repeating this method, we say that one family of integral curves (covering the total 
space) for the given set of parameters is equivalent to the other for another set of parameters, if 
there exists a diffeomorphism of space transforming the first family into the second. It is clear 
that there occurs the division into the classes of equivalence, and we are not very interested in 
differences between equivalent families. We argue that the corresponding changes in parameters 
do not alter anything on principle. So it is sufficient to demonstrate only one phase portrait for a 
given set of parameters underlining the features of the given class. 

However, for some critical relations between the parameters there occurs significant changes. 
These are the boundary relations of parameters, dividing, as usual, parameter space into regions 
of similar behavior. Thus accomplishes the qualitative classification of the mode of evolution of 
dynamic system. Now, giving the concrete value of parameters, we can define which region of 
parameters they correspond to, thus define the type of behavior. Moreover, given the specific 
initial conditions, we can answer the question to which region of phase space the evolution of the 
system lead in time. 
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In our cosmological model, numerical parameters A, a, B, /3 are related to the viscosity, while 
A and A are the (self)-coupling and cosmological constants. 

Initially, we consider the system of Einstein and Dirac equations. Solving these equations, 
we find the components of the spinor field and metric functions a, b, c in terms of volume scale 
T = abc of the BI universe. Finally, in order to find x from Einstein equations and Bianchi identity, 
we deduce three first order ordinary differential equations. Further for convenience we introduce 
a new function v inverse to X, i.e., V = l/x. 

The fact that the system has the dimension greater than 2, strongly complicates qualitative 
analysis. Note that well known Lorentz system of three ordinary differential equations with poly- 
nomial right hand side with degree less or equal to 2, possesses in some region of parameter space 
chaotic behavior known as a strange attractor and in that region there do not exist first integrals 
(i.e., globally defined invariants). Though the set of singularities is very simple, there exist only 
three singular (fixed) points: two focus and one saddle. The presence of such example does not 
allow us to make an optimistic conclusion on the basis of simple construction of our system (with 
polynomials in the right hand side and absence of singular points the in region of space we are 
interest in, which is even dynamically closed. 

Nevertheless, on the boundary of the the space £ = 0, as well as v = (t = +°°), which 
are dynamically closed themselves, the complete classification has been done. The dynamical 
closeness of these planes simultaneously as an obstacle for penetration from positive octant £ > 
A v > to the region with negative values. But, there are no singularities, fixed points (there are 
fixed points on the boundary) in the positive octant, we were not able to prove the simplicity of its 
behavior, e.g., presence of first integrals, as well as their absence. 

Thus let us go back to the system (2.1 1) in details. As it was already mentioned, tt is convenient 
to define a new function v = 1 jx. In this case the obvious singularity that occurs at x = vanishes 
and v = corresponds to T = °° while v = °° to T = 0. The system (2.1 1) on account of (2.13) 
takes the form: 

v = -3Hv, (3.1a) 
H = - — ft)) — (3H 2 - ke — A) + — (my + ^ - y _ n)2 ) , (3.1b) 

£ = 3H(3£H -co) +477 (3# 2 - K£ - A) -4t? [ K (mv + ^ + " )] ■ (3.1c) 
Let us now study the foregoing system of equations in details. 

A. Behavior of the solutions on v = plane 

As one can see, in this case the system (3.1) takes the form: 

H = |(3£#-<o) - (3// 2 -fC£-A), (3.2a) 
£ = 3H(3$H-Q)) +4r](3H 2 -K£-A). (3.2b) 

This system of equations completely coincides with the one when the BI universe is filled with 
viscous fluid only. The system in question was thoroughly studied in [18], hence we skip this 
study in the present report. 
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B. Behavior of the solutions on s = plane 



The plane e = is dynamic invariant, since £ | e=0 = 0. Depending on the sign of H this plane is 
either attractive or repulsive, namely, for H > it is attractive and for H < it is repulsive, since 

= -3H(l + C) <0. 
In presence of of spinor and scalar fields the system (3.1) at £ = has the form 



v = -3Hv, (3.3a) 
H = - 3 ^ + A + l( mv + ^^ ¥ ). (33b) 

The system (3.3) has the following integral curves 

Cv 2 - 

(Av' l + 1 



y2+n 

6H 2 = 2A + 2mV + Cv 2 -- (3.4a) 



where C is some arbitrary constant. 

The characteristic equation of nontrivial singular points on £ = plane for the system (2.11) 
takes the form 

2mA 2 v 2n+1 + 4AA 2 v 2 ' 1 + nv n+2 + 4mAv" +1 + 8AAv' 7 + 2mv + 4A = 0. (3.5) 



Depending on changes of signs in the sequence of A, m, A it has one, two or no solutions. 

In Tables Al, Bl, CI, Dl we illustrated the phase-portrait on £ = plane for a positive and 
a negative A, respectively for n = 1,2,3,4 and A < 0. In Tables A2, B2, C2, D2 we illustrated 
the phase-portrait on £ = plane for a positive and a negative A, respectively for n = 1,2,3,4 and 
A = 0. In Tables A3, B3, C3, D3 we illustrated the phase-portrait on £ = plane for a positive and 
a negative A, respectively for n = 1, 2, 3, 4 and A > 0. 

As it was mentioned earlier, here we deal with the multi-parametric system of ordinary nonlin- 
ear differential equation. In doing so we consider all possible variants independent to their physical 
validity. Therefore, we demonstrate the results obtained for a negative spinor mass (m < 0). 

The singular point around which the oscillation takes place has H = 0, and therefore, the trajec- 
tory of oscillation partially passes in the region which is attractive to the plane £ = and partially 
in the region that is repulsive. In the long run in the repulsive region at some moment the growth 
of £ becomes dominant. It results in the fact that £ becomes infinity within a finite range of time. 
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m < 



A<0 




A = 




A>0 



m = 




II 



m > 





g 



Table Al. Case with e = 0, n = 1 and X < 0. 



A<0 



A = 



A>0 



m<0 






m = 






m > 






g 



able El . Case with e = 0, n = 1 and A = 0. 



A<0 



A = 



A>0 



m<0 






m = 






m > 






Table A2. Case with e = 0, n = 1 



and A > 0. 
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m < 



A<0 



/, 




A = 



A>0 




m = 





II 



m > 





g 



Table Bl. Case with £ = 0, n = 2 and X < 0. 



A<0 



A = 



A>0 



m<0 





T 




m = 






m > 






g 



Table E2. Case with £ = 0, n = 2 and A = 0. 




Table B2. Case with £ = 0, n = 2 and A > 0. 
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m < 



/ 



A = 




A>0 




m = 





II 



m > 






g 



Table CI. Case with £ = 0, n = 3 and X < 0. 



A<0 



A = 



A>0 



m < 



m = 



m > 





g 



Table E3. Case with £ = 0, n = 3 and A = 0. 



A<0 



A = 



A>0 



1 — /" 



m < 




m = 





m > 





Table C2. Case with £ = 0, n = 3 and A > 0. 
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A<0 



/ , 




A = 




A>0 




m = 





II 



m > 




g 



Table Dl . Case with £ = 0, n = 4 and X < 0. 



A<0 



A = 



A>0 



m<0 



-r 



m = 



m > 



g 



Table E4. Case with £ = 0, n = 4 and A = 0. 



A<0 



A = 



A>0 



m < 




m = 





1/ 



m > 






Table D2. Case with £ = 0, n = 4 and A > 0. 
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C. Qualitative analysis of the complete system 

The system (3.1) in absence of viscosity, i.e., under r\ = and £, = possesses the following 
first integrals 



Fx 
F 2 



(6// 2 -2e-2A-2mv) 
V 2 



(3.6a) 
(3.6b) 



The second of them (3.6b) remains to be the first integral even after the introduction of bulk 
viscosity <§. The first one, i.e., Eq. (3.6a) under ^ ^ ceases to be the integral of motion. Never- 
theless, the introduction of bulk viscosity during the course of time generates definite displacement 
of the surface given by the formula (3.6a), which allows one qualitatively, i.e., based only on the 
continuity, compile the representation about the possible ways of evolution. 




5.0 7.5 10.0 





FIG. 1: Evolution of func- FIG. 2: Evolution of FIG 3: 3D view in 

tion inverse to volume volume scale v,H,e space 

scale 



v: 

1.5- 




FIG. 4: Evolution of func- FIG. 5: Evolution of FIG 6: 3D view in 

tion inverse to volume volume scale v,H,e space 

scale 



Harnessing the Tables 1, 2 and 3, helps one to understand the 3D phase portrait leaning on the 
continuous dependence of the velocity fields of the coordinates v,//,£ of phase space. 
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FIG. 7: Evolution of func- FIG. 8: Evolution of FIG. 9: 3D view in 

tion inverse to volume volume scale v,H,e space 

scale 




FIG. 10: Evolution of FIG. 11: Evolution FIG. 12: 3D view in 

function inverse to volume of volume scale v,H,e space 

scale 



In order to cover the infinite phase space completely, it is mapped on coordinate parallelepiped 
with its axes being the the arc-tangent of the corresponding coordinates. The lower horizontal 
plane always represents the e = plane. 

It should be noted that the introduction of spinor field notably complicates the evolution of the 
system. Contrary to the system in absence of the spinor field, the initial condition with H < 
already does not prevent in many cases thanks to the evolution of volume scale entering the half- 
space H > and thereupon, from the greater value of H repeats the evolution, approaching to the 
v = plane and displaying the classification from the table 1. In the vicinity of the borders £ = 
and v = the integral curves closely repeats the integral curves on the sides, each time at least to 
some extent. 

The general property of all the cases is the fact that in the half-space H > the velocity vectors 
are directed to the £ = plane, while in the other half opposite to it. As a result all he invariant 
curves fall on £ = 0, though not necessarily reach it. 

In the Figs. 10 - 3 we have illustrated functions inverse to the volume scale v(t) [Figs. 10,7,4,1], 
volume scale x(t) [Figs. 11,8,5,2] and phase portrait in v, H, £ space [Figs. 12,9,6,3], for a = 4, 
J3 = 1, £ = 1/2, A = 1, B = 1, n = 4, m = 4. 

Continuous and dot lines in the Figs. 1-12 corresponds to two different initial conditions. For 
A < depending on the sign of A there occur the following situations: (i) for X < there exist 
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separable plane which does not allow the solutions with initial condition in one part enter into the 
second one [cf. . 1,2,3, that corresponds to Table Dl-g]. (ii) For A > there is no separable plane 
in this case [cf. . 4,5,6, that corresponds to Table D2-g]. As one sees, a negative A, which is in 
fact the additional gravitational field, generates oscillatory regime of expansion. 

In case of A > there are only exponential regimes of expansion. For X < there is separable 
plane [cf. . 7,8,9, that corresponds to Table Dl-i], while for A > there is no separable plane [cf. 
. 10,1 1,12, that corresponds to Table D2-i]. 

IV. EVOLUTION WITH BLOW UP 

Studying the system of ODE let us imagine the integral curves in the space. It is very important 
to know the directional field given by this system. It is more important than the corresponding 
vector field. First of all, the integral curves, by definition, are tangent to the vector field, hence to 
the directional field, at those (peculiar) points, where vector field becomes trivial with the direction 
being indefinite. Secondly, like the vector field the directional field is also continuous (excluding 
the peculiar points), but it may be continuously continued at the boundary where the vector field 
might be infinity. 

We are interested in two aspects: how rapidly the solution can tend to infinity at the distant 
boundary (simply infinity) and how does it behave at the infinity. The way the problem is posed 
becomes reasonable when the space is closed by means of infinitely remote points in any given 
interpretation. 

We will follow the Pensele's ( ) principle of continuity - the properties of a system at continuous 
change from one common position to another without losing generality. We are interested in 
qualitative properties in solving the system of ODE. Deforming the vector field continuously, at 
the same time leaving the peculiar points unaltered, we don't change the qualitative behavior of 
the integral curves with an accuracy of topological equivalent. In this way we can simplify the 
analysis, substituting the initial system by a simpler one, constructed from convenient elements. 

A. Blow up 

The history of studying the regime with blow up is associated with S.P. Kurdiumov [23]. The 
study of the process of heat distribution in active and nonlinear medium led to an extremely dis- 
tinguishing feature, namely wave and localization. Mathematical models of demography detects 
critical moments: solution to the (time dependent) ODE may reach its limit within a finite time. 
The processes in the chromosphere of the sun possess a flashing (eruptive) character, but the mech- 
anism of energy transference does not detect the presence of predefined scale of time. 

To illustrate the detection of a characteristic time in the system with no explicit time- 
dependence, let us consider the following example. 

x=-x a ,xeR + (4.1) 

It has two solutions: a) x{t) = andb)x(t) = [jc(0) 1_a -t(l - a)]^. 

In case of b) the limiting value x{t*) = is reached at a finite time /* = \_ a , if oc < 1. Then 
both solutions mix up. At moment U the uniqueness condition (precisely, Lipshits condition) 
breaks down. 

The power law dependencies are typical for different types of catastrophes: from earth quakes 
and flood to stock exchange collapse and accidents in atomic power energy. 
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B. Infinity 

The joining of infinitely remote point to the space of ODE 

x = F(x),x£R + (4.2) 

we execute in the following way: let us make the change of variables x = |, s 2 + c 2 = 1. We call 
the point ±°° = A infinitely remote one. 

As a result we obtain a system of equations 

sc — cs = c 2 F fS 
ss + cc = 0, 

which on account of s 2 + c 2 = 1 leads to 



c- 



(4.3) 



s = c 3 i^), 
c = sc 2 F( s -) 



(4.4) 



Reducing the right hand side of the system to a common denominator in the vicinity of the point 
s = 1 , c = (but not on it) and then eliminating it, we do not alter the directional field. Preserving 
namely this meaning, we define the direction at this point. 

Let us go back to the system of equations and rewrite it in the form 

v = -3Hv, (4.5a) 
H = I(3£//-(e + p))-(3// 2 -£-A) + i<Mv), (4.5b) 
e = 3//(3^//-(e+ j p))+477 (3// 2 -e-A) -477<^ 2 (v), (4.5c) 

where 0i and 02 are the functions of T. 

In case of a spinor field only we have 0i (v) = mv + X(n — 2), fo.(v) = mv — Xv n . 

n+2 2 

Introduction of a scalar field gives 0i (v) = mv + 2 (i+Xv n ) 2 ' ^ 2 ( v ) = mV 2(i+Av") - 
Near the point £=«we make the following substitution £ = l//i. Then the system takes the 
form 

v = -3Hv, (4.6a) 
H = |5H/i-^ + i(l-C)/i _1 -3H 2 +A+^i(v), (4.6b) 

/i = 4A(-3H 2 +A + (j)2(v))Li 2 - a -9BH 2 li 2 -P +4ALi l ~ a +3H(1 + Qii. (4.6c) 

As it is seen from (4.6c) in the absence of viscosity (A = 0, B = 0) the blow up along the energy 
density is impossible. 

The answer, whether the blow up takes place in the past or in the future, depends on the sign of 
the coefficient at jl with the lowest power. 

Let A = 0. In order to the blow up takes place at finite H, it is necessary that /3 > 1 . In this case 
the singularity will be in the future, i.e., we have Big Rip. 

Now consider the case with B = 0. In this case the blow up takes place in the past (Big Bang) 
if a > 1 . 

In the figures illustrated below we plot the trajectories on which the infinite energy density e is 
achieved in a finite time. The blue line indicates past while the red one the future. 

In the figures 14 and 13 we show the evolution of T, H and £ relative to each other. In both 
cases there exists possibility for infinite growth of energy density at infinitely large volume, i.e., 
there occurs so-called Big Rip. 
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FIG. 13: The trajectory of evolution in case of an interacting 
spinor and scalar fields with a = 4,j6 = l,£ = 1/2, A = l,B = 
1 ,m = 4, A = -1,A = -1 



V. CONCLUSION 

Recently a self consistent system of nonlinear spinor and gravitational fields in the framework 
of Bianchi type-I cosmological model filled with viscous fluid was considered by one of the au- 
thors [19, 20]. The spinor filed nonlinearity is taken to be some power law of the invariants of 
bilinear spinor forms, namely I = S 2 = (\j/\j/) 2 and J = P 2 = (ixj/y 5 \]/) 2 . Solutions to the corre- 
sponding equations are given in terms of the volume scale of the BI space-time, i.e., in terms of 
T = abc, with a, b, c being the metric functions. This study generates a multi-parametric system of 
ordinary differential equations [19, 20]. Given the richness of the system of equations in this paper 
a qualitative analysis of the system in question has been thoroughly carried out. A complete qual- 
itative classification of the mode of evolution of the universe given by the corresponding dynamic 
system has been illustrated. In doing so we have considered all possible values of problem param- 
eters independent to their physical validity and graphically presented the most distinguishable in 
our view results. 

The system is studied from the view point of blow up. It has been shown that in absence of 
viscosity the blow up does not occur. It should be emphasized that phenomena similar to one in 
question can be observed in other discipline of physics and present enormous interest from the 



Anisotropic cosmological models with spinor and scalar Gelds 




FIG. 14: The trajectory of evolution in case of a spinor field with 
self-action at a = 4,/3 = 1,£ = 1/2, A = 1,5 = l,m = 4, A = 
-1,A = -1 



point of catastrophe, demography etc. 
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